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Wc propose theoretically a method for the control of the frequency of the nonlinear optical 

response of a model nonlinear medium driven by an electric field and suggest an experimental 

realization on the basis of isotropic LiF, or NaCl crystals and pulsed Nd glass laser. The 

theoretical background of the proposed optical control is the nonlinear resonance and the 

Cl^ resonance overlapping which is typical for nonlinear systems. It is shown that by a proper 

O 

choice of the parameters of the applied pulses a diffusive growth of the oscillation amplitude 
can be achieved resulting in a controlled switching of the system's frequency. The maximum 
>^ frequency of the nonlinear response is identified to be proportional to the intensity of the 

applied noise. An experimental suggestion for particular materials is proposed. 

PACS numbers: ... 



I. INTRODUCTION 

With the accumulated volume of information that has to be transmitted, assorted and stored 
swiftly and efficiently there is growing demand to utilize for these purposes methods and techniques 



from nonlinear [lHJ] and fiber optics jo-tlC| 



For an optical system that is slightly off its equilibrium state, the model of a linear oscillator 
is a well-established picture that yields a good qualitative and quantitative agreement between 
theory and experiment. In particular, for different types of optical media, optical characteristics 
such as the electric susceptibility, the polarization and the transmission-refraction coefficients were 
predicted correctly by the linear oscillator model [l-4]. However, systems that are far off the 
equilibrium possess usually physical features that cannot be explained using linear models. For 
example, the phenomenon of the self-induced transparency and the generation of higher harmonics 
were successfully explained by a nonlinear oscillator model [fjj]. The problem of the light-induced 
critical behavior in nonlinear systems is a classical problem and was addressed in many studies 



2 



llUl3l]. An interesting phenomenon observed is the polarization related optical bistability that 
occurs due to the dependence of the transmission and the absorption coefficients on the oscillation 
amplitude that may well be chaotic for a nonlinear system. The polarization dynamics is usually 
defined by the Kitano map fl3 ] and undergoes a bifurcation transition with the increase of light 
intensity. However, the problem of the optical transistor we are going to address in the present 
study is not connected with the polarization dynamics, but with the controlled switching of the 
frequency. The physical mechanism for the functioning of the proposed optical transistor is the 
phenomenon of nonlinear resonance and resonance overlapping. This phenomenon is typical for 
externally driven nonlinear systems and can be utilized for the diffusive control of the frequency 
of system. For nonlinear optical media we may expect a response with the frequency being several 
times higher than the frequency of the applied external driving field. The question that arises 
naturally is whether it is possible to tune and control the outgoing frequency. Here, we develop 
an efficient frequency switching scheme which can be easily realized experimentally. We focus 
on reaching the maximum desired outgoing frequency that defines the maximum transmission 
bandwidth for the optical transistor. We also address an even more important issue - the freezing 
of the target frequency fixed which naturally required for the stability of the device. 
The paper is organized as follows: In the second section we specify the problem for the case of 
a single oscillator model and evaluate the maximum possible response frequency which can be 
obtained from the optical media using a monochromatic external driving field. We also consider 
the scheme based on the polychromatic pumping technique. We utilize the method of nonlinear 
resonance in order to derive the corresponding Fokker-Planck equation that describes the evolution 
of the system frequency when the oscillator dynamics is chaotic in time. Additionally, we study the 
consequences of both dissipation and thermal noise. In the third section we generalize the model 
for the case of interacting oscillators paying a special attention to the collective dynamics and the 
correlation effects between different oscillators. 

II. SINGLE OSCILLATOR MODEL 
A. Model 



Generally the oscillator model views the active medium as consisting of a set of classical an- 
harmonic oscillators in a unit volume. Each oscillator stands physically for an electron bound to a 
core or an infrared-active molecular vibration. The total energy of the electron in the presence of 



a driving electric field reads 



in 



H = H {x,p) + H nl (x) + eV(x,t), 

'0 



tfo = f + ^- (1) 
2m 2 

i?NL = —ax 3 — /3a; . 

Here (^p) is the energy of the harmonic oscillator, .Hnl is a nonlinear correction of energy, and 
eV(x, t) is the interaction with the electromagnetic field. In particular we will consider two possible 
cases for the external driving: a) an applied external monochromatic electric field eV (x, t) = 
xVo(t), Vo(t) = e/cos(Jli), with the amplitude / and frequency fi, and b) an infinite series of 
applied delta pulses V (x,t) = xVq J2^=-oo $(t ~~ n T), Vq = ef, with the amplitude / and the 
interval between the pulses being equal to T. In Eq. ([1]) e and m are the charge and the mass 
of the electron, a and j3 are the nonlinearity coefficients, x and p are the coordinate and the 
momentum of the electron, and ujq is the eigenfrequency of harmonic oscillations. In what follows 
we assume that the eigenfrequency ujq and the driving frequency of the monochromatic field Q, are 
optical frequencies, while the period between the pulses is chosen to be of the order of a picosecond. 
Therefore, the following condition is met <C T. In what follows, we set time interval between 

pulses T as the unit of the time and therefore T = 1. These parameters are typical for the Nd- 



gas laser [j] operating in the mode-locking regime. The mode 
with a nonlinear oscillator model seems a great simplification 



ing of the driven physical system 
111. However, as mentioned in the 



introduction it has been established as a very efficient and to a large extent reliable approach in 

Q. 



nonlinear optics 



B. Monochromatic driving 

In the case of a monochromatic driving, we add to the equation of motion (JTJ) a phenomenological 

decay term 72; responsible for thermal loses in the system (pQ) 

2 3a 2 4/3 3 7 

x = — UnX x x x H cos(sZr). (2) 

m m m m 

In the case of a weak nonlinearity corresponding to the state that deviates slightly from the equilib- 
rium, i.e. (xU)) max < ^y^F, equation (2) can be solved analytically using canonical perturbation 
theory 



161 ] . Consequently, for the amplitude frequency characteristics one obtains the following 



expressions 



/ 3(3 5a 2 
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FIG. 1: Schematically shown dependence of the oscillation amplitude A max [15| on the detuning between 
the frequency of the external field and the frequency of the linear oscillations Alo = Q — ujq. The maximum 
value of the amplitude belongs to the unstable domain BC of the hysteresis loop. 



In case of weak nonlinearity maximal values of the frequency shift corresponds to the case a = 0: 

3/3 



Ulr, 



^0 + 



— A 2 

8u m max ' 



(3) 



where A max is the maximum square of the amplitude which may be achieved during the oscillations 
for the case a = 



V 2 
Auh 2 ' 



(4) 



Since the maximum values of the oscillation amplitude is limited by the condition @, the maximum 
nonlinear frequency shift 3/3/(8cJo"i)^4 max an d> consequently, the maximum frequency w max that 
can be achieved using the monochromatic driving external field is limited too. Another drawback 
of the monochromatic driving is the fact that the maximum frequency w max belongs to the unstable 
domain of the amplitude frequency characteristics (cf. Fig. [1]). Hence, even a small fluctuation of 
the driving frequency Aw = Q — ujq leads to a drop of the maximal frequency w max down to the 
linear frequency loq. 



C. Dynamical stochasticity and diffusive control of the frequency 

Looking for a more efficient frequency switching scheme for the optical transistor, we con- 
sider an infinite series of delta-function pulses (kicks) applied to the system, i.e. V (x, t) = 
x VqT Y^n=-oo ~ n T), where Vq = ef. First, we neglect dissipation in the system. Transforming 
to the canonical variables (I, 9) of the harmonic oscillator (with frequency ojq), i.e., x = J cos 9 



and p = y/2IujQm sin 9, we deduce from eq. (1) that after averaging over the fast phases [ijj], the 
relations apply 

H = H (I) + V(I,9,t), 

H (I)=I<j + torp(—) , 
\muj j 

oo 

V(I,9,t) = V(I,9)T £ 5(t-Tn), (5) 



where V(J, 0) = Vb(J) cos 6* and V = VfoCO ~ Vb) 

67T/3 

£j(i) = w H 5 2 

1 11 '"^0 
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From the system of eqs. Q we see that the frequency of the system depends on the action oj{I). 
Consequently, when controlling I, we can control the frequency of the system. The Hamilton 
equations of motion corresponding to the energy given by eqs. (|5]) read 

f _ c dv(i,e,t) 

." " d9 dv(iet) ^ 

ol 

The benefit of using the short kicks is the simplicity of the corresponding equations of motion 
([6]). Since the motion is governed by either kicks or by the free evolution inbetween two kicks, the 
equations of motion ([6]) can be integrated straightforwardly over one temporal period by splitting 
the evolution operator in two parts T = TrTs- Here 

f R (I,9) = (I,9 + oj(I)T) (7) 

is the evolution operator responsible for the free rotations between the pulses, while the operator 
T$ describes the influence of the pulses 

ft{/l9)= ( / . eT TO « + «T^). (8) 

Taking into account eqs. (|7|) and ([8]) and assuming that Qit f3 / (m 2 ujq) 3> Vn I (2J mujo) , the solution 



of the equations of motion ([6]) can be found in the form of the standard map 17] 

I n +i = I n + K sin 9, 



(9) 

n+l — UIqT + 9 n + I n +l, 



where K a V^T 2 . 



Despite its relative simplicity, the standard map displays a complicated and rich behavior 

0- 

For large K > 1 the phase space of the system is predominantly chaotic and is an invariant 



Kolmogorov-Arnold-Moser (KAM) torus [171 . Il8l ] . In the chaotic sea that covers most of the phase 
space, the dynamics of the action variable I is expected to be diffusive in time. Taking into 
account the following values of parameters that corresponds to the Iceland spark (Iceland crystal) 
and the isotropic crystals LiF, NaCl and Nd glass lasers T = 10" 12 [s], ujq = 10 15 [s- 1 ], 

m = 10~ 31 [kg], e = 10~ 19 [C], /3 = 4 • 10 11 [kg/(m 2 s 2 ], x max (0) ~ 1(T 8 [m], we see that the 
condition K > 4 imposes the following restriction on the driving field amplitude 

With increasing K the principal resonance island gradually becomes smaller and disappears for 



about K > 5 



16]. When the condition K > 1 holds, then the dynamics of the action variable is 

n 

chaotic and is described by the following Fokker-Planck equation [17|] 

where A = (A/), D = {^(A/) 2 ^, AI = I n +i — I n - Taking into account eqs. ([9]) and (fTT|) we find 

,4 = 0, D = TUVsmey) = T— . (12) 



The Fokker-Planck equation can be simplified as 

dF _ 

and finally we obtain 



dF 1 d 2 F . . 

-D^rl (13) 



(P) = v// 2 + Dt. (14) 

Using the connection between the action variable and the nonlinear frequency © and (|14p we 
derive an explicit expression for the time dependence of the nonlinear frequency 

<u;(I)} = U (y/{P)) =u + JLy/q + Dt = co + ^I VT+Wi. (15) 

Here D' = ^572-, = UJ o mx ma x ~ 10~ 32 [Jsl. Taking the values of the parameters from Ref. [ul ]. 
i.e. T = 10" 12 [s], u = 10 15 [s-1], m = 10" 31 [kg], e = 1(T 19 [C], f3 = 4 • 10 11 ]kg/(m 2 s 2 )], 
^max(O) ~ 10~ 8 [m], we arrive at the following estimate ^i^/o ~ 24vrw , D' ps 10 12 [s^ 1 ] = 1/T. 
The meaning of the expression (|15p is the following. Imagine, that the pulses of a light beam are in- 
cident on a nonlinear optical medium which is modeled as an ensemble of nonlinear non-interacting 
oscillators. The light heats up the oscillator motion and converts the regular small oscillations with 
the harmonic frequency ujq into the irregular oscillations with the nonlinear frequency An 
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averaging using the Fokker-Planck equation is equivalent to the averaging over an ensemble of 
non-interacting chaotic oscillators. As a result, the averaged quantity < uj(I) > exactly matches 
the frequency of the resulting outgoing signal that could be expected from optical media in an 
experiment for real physical materials. For example, for Iceland spark and isotropic crystals such 
as LiF nad NaCl, the model of a driven nonlinear oscillator gives the theoretical estimates that 
are in a good qualitative and quantitative agreement with the experimentally observed physical 
properties 0, 4]. Suppose, that our aim is to achieve a signal from an optical medium with the 
frequency that is a multiple of the initial harmonic frequency ujq, Nujq, where N is an integer 
number. According to Eq. (|15p . as a consequence of the stochastic heating, the nonlinear optical 
medium responds at uj = Nujq after a time defined via the following equation 

1 



(16) 



Eq. (|16j) corresponds to the normal diffusion case. In the case of a super diffusion it modifies to 



1 



, l<n<2 (17) 



The analytical result given by eq. (|16|) is in a good agreement with the exact numerical calculations 
for the map given by eq. (j2J) (cf. Fig. [2]). Inserting standard values of the parameters (listed below 
the Eq. (22)) in the Eq. (23) provides a very simple estimation for the frequency switching time: 
t c = 10 4 r(^±) 2 , < u(t) >= ujq + uolQ- 2 Ky/l + D't, >/<■*"(*)> =y /l + D't, D' = 1/T. 



Now we immediately can estimate time t c (number of kicks, since unit of time is the interval 
between pulses) for the switching of initial frequency loq to the final frequency uj = Nujq, for a 
given stochastic parameter K. 

Indeed, the numerical calculations confirm the diffusive increase of the action variable for this 
case. 

Note that in eqs. (|16p and (|1T|) we do not set any limitations on N. This means, that one can 
reach a final frequency as large as one wishes in principle. This is a principle benefit of the applied 
pulses, since in the case of a monochromatic pumping we have a maximum frequency limit (eq. 
([3])). Nevertheless, the disadvantage of eqs. (|16p and (|17p is the frequency stabilization problem. 
Due to the diffusion (by the nonlinear nature of the problem) the frequency will increase even after 
reaching the desired value uj = Nujq. However, as will be shown below, for more realistic models 
which include the decaying and the thermal effects we can circumvent this problem. 
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FIG. 2: /„ dynamics in the dimensionless units ^ <I j^ > - = Vl + D't = y/1 + t/T on the scale n <G [0 : 1000] 
for 7 = 0. a) and b) correspond to the stochasticity coefficient K = 0.5; c) and d) describe the dynamics for 
K = 5.0. < 1% > x / 2 (b, d) is defined as the ensemble average over 1000 realizations of the initial random 
values from the intervals I n =o €(0:1) and 6 n =o <= (0 : 2n). 

D. Dissipative map and frequency freeze 



Obviously, a realistic model should include not only the anharmonicity of the oscillations but 
also the dissipation that may happen in optical media due to thermal losses. We therefore generalize 
our model (eq. by adding a phenomenological damping term —-yx. Following the standard 

procedure given by eqs. (J7J) and (JHJ) and splitting the evolution operator in two parts T = TpT/j, it 



can be shown that the dissipation term modifies the standard map in the following way 

I n+ i = e"°7 n + if (7) sin6> n , 

$n+l — #n + In+U 



17] 



(18) 



Here I = (I — and 7 = 7T are the rescaled dimensionless action variable and the rescaled 

decay constant, respectively. The Fokker-Planck equation has the form of eq. (jlip . however, the 
coefficients given by eq. (|12p take on a different form 



.4 



(1-e" 7 )/, D= (1-e-^) 2 I 2 + ^K 2 e~^ 



(19) 
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Taking into account eq. (|19|) . the Fokker-Planck equation reads 



« -( 1 -^'^< /F ™ + 2^U 2F < / »)J + 4 /f e (20) 

Eq. (|2U|) is quite involved and therefore finding of exact analytic solution for distribution 
function is complicated problem. However for evaluation of the second moment of random variable 
< I 2 >= Jq 00 dl • I 2 F(I,t) we do not need explicit solution for distribution function. Multiplying 
both sides of eq. (I20p on I 2 , integrating over I and taking into account f£° 9F qP I 2 dl = < I 2 > 
we deduce: 

I 2 ) = I 2 exp (-2 7 i) + ±-K 2 [1 - exp (-2 7 t)] , (21) 
with the asymptotic limit 



/-/ \ 2 K 



2^7 

Consequently, for the nonlinear frequency we obtain 



V \ J o / 2^/7 (22) 

(A/2) hK 



6irP K 2 K 2 

Wmax (t ->■ oo) = w H 5— o W 1 + — ~ ^o + 247rw W 1 + — . (23) 

m^Wg V ^7 y 47 

From eq. (|23|) we infer that in the case of dissipative systems, the maximum frequency w max (i — > oo) 
has saturated values that can be controlled by tuning the chaos parameter K. The analytical results 
are in a very good quantitative agreement with the results of exact numerical calculations (cf. Fig. 

ED- 

E. Noise-induced enhancement 

To finally close the consideration of the single oscillator problem, we discuss the effect of the 
controlled noise applied to the system. We consider a modification of the interaction, when the 
fixed kick amplitude is replaced by a step-dependent random amplitude and instead of eqs. (jT8j) 
we obtain now: 

l n +\{m) = e~7I n + if (7, £ re ) sin6> n , 

On+i(m) = 9 n + J n +i, (24) 
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FIG. 3: /„ dynamics on the scale n £ [0 : 1000] for K = 5.0 and various 7. v/< 1% > is denned as the 

ensemble average over 1000 realizations of the initial random values from the intervals Iq £ (0 : 1) and 

X 

O £ (0 : 2tt). For finite 7 the theory predicts V '< ^ > = ^ I*-kx> (eq. [22]). For 7 = 0.005 and 7 = 0.01 

2 V7 



the analytical expression yields y/ < 1% > ~ 25 and yj < 1% > ~ 35, respectively. The numerical calculations 
illustrated by the figure for 7 = 0.005 and 7 = 0.01 visibly saturate to those values. The time needed to 
approach the saturation value scales as « I/7. 



where £ n is a Gaussian white noise = 0, with the noise strength q (£(£ + At)£ffl) 



q5(At) which is implemented into the numerical procedure via standard methods, see e.g. 



2m- 



It is well-known that noise is sometimes profitable and can indeed enhance signals in many 



cases 



211 ] . We intend to explore the consequences of the applied noise on the transistor frequency. 
With this purpose we numerically integrate the recurrence relations (|24p for different strengths of 
the applied noise. The results of the numerical calculations are presented in Figs. |4] and [5j 
From Fig. [5] we see that the saturation time is noise-independent while the saturated values of 
action increase with the noise strength. An increase of the frequency upon enhancing the noise 
is also obtained in other physical phenomena. Particularly, in experiments using the ferromag- 



netic resonance (FMR) technique 



221 ] the resonance magnetic field, which is associated with the 



resonance frequency, grows with the increasing noise level. Those experimental findings are well 
reproduced by numerical simulations including thermal noise 2314251] . 



In addition, the effective stochastic coefficient K turns out to be proportional to the noise strength 
q. This is obvious from the relation (|22p y/ (I 2 ) = and the numerical data presented in Fig. 
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[5j We obtain a linear function of the noise K(q) = K + ctoQ, where the constant ao can be ex- 
tracted from Fig. [5j Consequently, for the maximum frequency of the optical transistor we have 
the following estimation 



^max (t -> OO) « W + 2AtTU J 1 + - 



47 (25) 



Now we have all ingredients concerning the optical transistor factor. It looks relatively simple and 
should easily be accessible in the experiment 



N = = ! + 24 Jl + (^ + "Dg) 8 , (26) 

w V 4 7 

The only control parameters are the parameters of the applied pulses: the time interval between 

the pulses T, the amplitude of the pulses / and the intensity of the applied noise. 

Before generalize our approach for the ensemble of oscillators, we will consider role of fluctuations 

in order to have clear understanding how precisely switching procedure can be performed. Note that 

nonlinear frequency is a linear function of action u(I) (5). Along with the diffusive growth of mean 

value of action (I) actual values of action fluctuates, leading to the fluctuation of the frequency oj(I). 

Maximal values of the deviation of action was evaluated in and reads (5I) max » y/v^y 1 . 

Comparing maximal values of fluctuation (SI) with the maximal values of the action increment 

observed for the dispersive map (22) we deduce ^ S ^ max Typical values of the dimensionless 

dispassion constant is small 7 << 1, while stochastic parameter is large K > 1. Since nonlinear 
frequency is linear function of action for the frequency we have the same estimation ^jsai 

'Alu) 2 



( 5 ^)max — 1 Here (Sui) max is the maximal fluctuation of the frequency due to the 

fluctuation of action, and J (Au>) 2 is the maximal diffusive growth of the frequency. Consequently 
we may conclude that effect of fluctuations is reasonable small and accuracy of the frequency 
switching is high. 



III. ENSEMBLE OF OSCILLATORS 

Here we generalize our approach to a chain of interacting oscillators. In non-integrable Hamil- 
tonian systems with a large number of degrees of freedom F, the invariant KAM torus cannot 
isolate the stochastic layers formed in the vicinity of the separatrix lines 



2614311]. Naively, one can 



assume that the large number of degrees of freedom leads to stronger chaos. In principal, such a 
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FIG. 4: /„ dynamics on the scale n <E [0 : 1000] for fixed K = 5.0 and 7 = 0.005. y/ < 1% > is denned as 
the ensemble average over 1000 realizations of the initial random values from the intervals /q 6 (0 : 1) and 
6*o € (0 : 27r). The curve with q = 0.0 corresponds to that shown in Fig. [3] for the same 7. q is the noise 
strength, i.e. + At)£(f)) = q(5(At). 
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FIG. 5: y/ < 1^ >|t->oo(<?)-dependence obtained as a linear fit for regions where \J < 1% > saturates. The 
approximate values of a for the corresponding increasing damping parameter are found using eq. (|25|) and 
are a (7 = 0.001) w 19, a (j = 0.005) « 9 and a (7 = 0.01) w 7. 
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statement is correct, however, even in the presence of strong chaos the phase space of the system 
still may contain regular islands. However, as demonstrated in Ref. 28J the measure of the regular 



islands decays exponentially with the dimension of the system. For real materials, we assume that 
F>1 and consider the effect of dissipation and noise as in the F = 1 case in Sec. II. Keeping in 
mind the experimental situation, we will pay a special attention and precisely evaluate the time 
dependence of the diffusion coefficient. We assume that the interaction between the oscillators is 
mediated by an external driving field and generalize expression ([5]) for the multidimensional case 

N 

H = J2(H o (Ii) + V(I i ,0 i ,t)), 

i=l 

I \ 2 
1 \ 



H (I i )=I i u + torp — i- , 
\muJoJ 

oo 

V(I u 6 i ,t) = V(I i ,6 i )T J2 S(t-Tn), (27) 



where V(I U 9$ = V (Ii) cosft, V = V Q {h) « V 0] 

i t \ 6-7T/3 

U)(Ii) = UJq H 5— o H- 

m z ujQ 

Eqs. (|27p generates the map that can be considered as a generalization of the K. Kaneko, and T. 
Konishi model [2j 



P n+1 = r n + K S m6\-a- sin (o^ 1 - + a • sin - ^ , 
= uoT + 9 l n + r n+1 . 



(28) 



In the linear limit when the deflection of the angles from the equilibrium position 9 l n is small 
and sin#^ can be linearized, the system (|28|) has a solution of the type of localized breathers 



331 ] . However, we are interested in the strongly nonlinear chaotic case. The results of a numerical 
integration of the map given by eq. (|28p are presented in Fig. [6l which shows that, when increasing 
the coupling strength between the oscillators, the diffusion process is indeed enhanced. 
Following the standard approach already used above and including dissipation we obtain the new 
map 



4+i = e" 7 /; + K(i) sh< - a( 7 ) ■ (W(^ +1 - 9\ 

ai+l ai _j_ ri 

*M = i^f) K - 

a(7) = f 1 ^ J 0(7)1 



(29) 



where 0(7) is the oscillators interaction strength. 

Results of numerical integrations of the dissipative map given by eq. ([29]) are presented in Fig. [71 
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FIG. 6: / n -dynamics of a large chain of 1000 oscillators on the scale n 6 [0 : 1000] for K = 5.0, 7 = 0.0, 
£ = 0.0 and various interaction strengths a. The initial configuration is randomly chosen from the intervals 
Jj £ (0 : 1) and 8 l Q <G (0 : 27r). T„ is defined according to Table Q] 

The time behavior of the diffusion process for the case of coupled oscillators is qualitatively the 
same as for the case of a single oscillator. Namely, the saturation of the diffusion is observed. The 
novelty is, however, that the saturated value of the action grows with increasing coupling between 
the oscillators. 

By analogy to the standard map of eq. (124D . the thermal noise can be included as follows 



Numerical results from Eqs. (|3Q|) are presented in Fig. [HJ 

We see that in the case of ensemble of oscillators thermal noise again results in increased 
saturated values of action. However absolute saturated values observed for the ensemble of 
oscillators are bit larger than values corresponding to the single oscillator case see Fig. 
Therefore we can conclude that nonzero interaction strength a(j) (cf. Fig. [7J enhances the 
diffusion process. 



r 



n+l ~ 



e-rP n + K{i) sin^ - a( 7 ) ■ (sin^ 1 " <£) " sin(^ - 9^ 1 )) 



n+l' 




(30) 



15 




200 



400 600 
Time step n 



800 



1000 



FIG. 7: ^-dynamics of a chain of 1000 oscillators on the scale n E [0 : 1000] for K = 5.0, 7 = 0.005, £ = 0.0 
and various interaction strengths a. Initial configurations are as in Fig. [6] /„ is defined according to Table 
ffl 



We also numerically evaluate diffusion coefficient using the following expression 

N 



D n = D(t = nT) 



1 N 1 

-T-\r 



(31) 



When the diffusion in the phase space is normal, there exists a finite constant = lim^oo D(t). 
Otherwise, in the case of an anomalous diffusion we have a polynomial decay = l/t s , 5 > 1 

El- 

The results of the numerical integration for eq. (|3ip are presented in Fig. [9j We infer from the 
figure that in the absence of the dissipation the diffusion is normal and the diffusion coefficient is 
a constant. The presence of the dissipation leads to a decay of the diffusion coefficient, confirming 
thus the saturation of the outgoing frequency. Hence, the considered chain of nonlinear oscillators 
allows for a controlled energy absorption, and may act as a prototype of a transistor in the optical 
frequency domain. 



IV. CONCLUSIONS 



The main finding of the present study is that the maximum frequency of the optical transistor 



is given by expression (|25j) . i.e. w max (t — > 00) ujq + 2A7ru!o\ 1 + 



(K + apgY 
47 



where cjq is the 
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FIG. 9: Time propagation of the diffusion coefficient D n given by eq. (|31[) in the non-dissipative and 
noise-free case (thick black curve), in the case of small dissipation 7 = 0.005 (thick red curve) and in the 
case where both dissipation 7 = 0.005 and thermal noise q = 1.0 are present (thin blue curve). The inset 
demonstrates the effect of growing thermal noise. 
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initial linear frequency, and the coefficient K is the coefficient of the stochastisity, which can be 
easily controlled by the amplitude of the applied pulses / and the intervals between them T. The 
coefficient 7 describes the phenomenological decay constant, q defines the strength of the noise and 
cko is a numerical factor which decays with increasing dissipation in the system ao(7 = 0.001) ~ 19, 
^0(7 = 0.005) « 9 and «o(7 = 0.01) « 7 (See Fig. 6). The observed switching time is equal to 
about 100T, i.e. T ~ 10~ 12 [s]. The amplitude of the applied pulses should satisfy the condition 
/ > 10 7 [V/m], which is realistic for the Nd glass lasers 0, [J. 

In addition to a single oscillator model, we have also considered a model of coupled oscillators, 
where for both models the inclusion of noise into the applied electric pulses is shown to be an 
efficient tool for increasing the maximum frequency of the optical transistor. 
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Appendix 



TABLE I: Definitions employed for calculations. 



Single body problem TV = 1 


Many body problem i = 


= 1,...,N, N = 1000 


• one realization: 

J n — V 1 n 


• one realization: 

»=1 i=l 


• £ 7^ 0, many realizations: 


• £ 7^ 0, many realizations: 


V< In > = \ 


1 R 9 

r=l 


1 N 


1 R 9 

r=l 


R = 


1000 


R = 1000 
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